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The duality and other symmetry properties of the effective conductivity a e 
of the classical two-dimensional isotropic randomly inhomogeneous heterophase 
systems at arbitrary number of phases N are discussed. A new approach for 
a obtaining <r e based on a duality relation is proposed. The exact values of 
<T e at some special sets of the partial parameters are found. The explicit ba- 
sic solutions of the duality relation, connected with the higher moments and 
satisfying all necessary requirements, are found at arbitrary values of partial 
parameters. It is shown that one of them can describe <r e for systems with a 
finite maximal characteristic scale of the inhomogeneities in a wide range of pa- 
rameters. The other solution, connected with a mean conductivity describes <x e 
of the introduced random parquet model of N -phase randomly inhomogeneous 
medium in some mean field like approximation. The comparison with the known 
effective medium approximation and crossover to the continuous smoothly inho- 
mogeneous case are also discussed. 



PACS: 72.80.Ng, 72.80.Tm, 73.61.-r 

1. Introduction 

The properties of the electrical transport of classical macroscopically inho- 
mogeneous (randomly or regularly) heterophase systems, consisting of N (N > 
2) phases with different conductivities Oj (i = 1,2, ...,N), always have had a 
great significance for practice (see, for example P3G])- Their importance revived 
in the last decades under investigation of various new materials, but now the 
inhomogeneities have more smaller scales (from microscopic and mesoscopic till 
macroscopic) E] ■ The main problem in a theory of the electrical transport 
of heterophase systems is a calculation of the effective conductivity a e at ar- 
bitrary partial conductivities and phase concentrations. Though some general 
results have been obtained for weakly inhomogeneous media PP, this problem 
remains unsolved for inhomogeneous media with strong inhomogeneous fluctu- 
ations, which can be connected with strong irregularities as well as with large 
differences of partial conductivities, due to a difficulty of the averaging over 
such fluctuations. In order to obtain any explicit expression for tr e a mean 
field like procedure, the effective medium (EM) approximation, has been of- 
fered many years ago |S] . It is applicable for all dimensions and describes well 
er e of heterophase systems for weak and moderate inhomogeneities for systems, 
admitting the random resistor network (RRN) representation [HJ . Unfortu- 
nately, the RRN model, being a discretised version of the continuous medium, 
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is not a good representation for all continuous media with moderate inhomo- 
geneities Moreover, already in 3-phase case, the EM approximation gives 
for er e a complicated algebraic expression, which is practically untractable (see 
below section 4 and Appendix C). For this reason one needs other methods 
and approximations, which will be more adequate for continuous random inho- 
mogeneous media. Fortunately, in two dimensions there is another possibility 
to treat the problem. As is known, the two-dimensional systems have an exact 
dual symmetry, which is a consequence of a duality between potential and diver- 
genceless fields in two dimensions. In case of two-phase isotropic self-reciprocal 
systems it gives the exact relation, connecting a e at the inversed values of the 
partial conductivities (Tj, [TrjlITT] 

a e (ai,a 2 ,x)a e (al/ai,aQ/a 2 ,x) = (Tq, (1) 

where Co is some arbitrary parameter with a dimensionality of a conductivity, 
a; is a concentration of one (for example, the first one) phase. For self-dual 
systems (i.e. systems, satisfying (1) and symmetrical under the permutation of 
the partial parameters (a±,x±) <-> (02, £2) ) the relation (1) can be represented 
also in the form, connecting a e at adjoint concentrations x and 1 — x 

a e (ai,a 2 ,x)a e (a 1 ,a 2 , 1 - x) = aicr 2 . (2) 

This duality relation (DR) allows to find the exact universal, remarkably simple, 
formula for a e at equal phase concentrations x\ = x 2 = x c = 1/2 

<J e = \/o\o 2 - (3) 

Recently it was shown that the DR with some additional assumptions can be 
used as a basic relation for a determination of er e of inhomogeneous two-phase 
systems, and two different explicit approximate expressions for a e have been 
derived [51ll2j. They differ from the EM approximation and satisfy all necessary 
requirements. One of these expressions, generalizing the exact Keller - Dykhne 
formula (3) on arbitrary concentrations 

<T e (0-1,02,2;) = afcr^ x , (4) 

describes the effective conductivity of self-dual systems with compact inclusions 
of one phase into another and is applicable almost everywhere, except the re- 
gions, where one of the partial conductivities <7j — > The other expres- 
sion for cr e describes a e of two-phase random parquet model with layered (or 
striped) structure. The existence of various expressions for a e means also that 
it is a nonuniversal function, depending on the structure of the inhomogeneities. 

In this paper we will study general properties of a e of two-dimensional self- 
dual heterophase systems with arbitrary number N of phases and show that 
the analogous expressions can be obtained for systems with N > 2 . Firstly, 
we will discuss general properties of these systems, including a duality rela- 
tion (a section 2), and find the fixed points of the duality transformation and 
the corresponding exact values for a e (a section 3). Then, in the section 4, 
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we consider the EM approximation for N -phase systems and show how the 
found exact values can be reproduced in this approximation. In the section 5 
we propose a general solution of the duality (or inversion) relation and obtain 
some explicit expressions for er e , two of which generalize the abovementioned 
solutions of N = 2 case. Then, two physical models of continuous random 
inhomogeneous media, whose effective conductivities coincide with these two 
expressions, will be presented in the sections 6 and 7. The first model repre- 
sents the N -phase systems with compact inclusions of different phases. For 
this model we deduce in the framework of the finite maximal scale averaging 
approximation (FMSA approximation) introduced in previous papers I12| 
the functional self-consistency equation for the effective conductivity and find 
a solution of this equation, its physical meaning and properties. The second 
model, generalizing the two-phase random parquet model on N -phase case, 
gives in a mean field like approximation another expression for effective conduc- 
tivity, coinciding with the second solution of the duality relation. These results 
demonstrate that the effective conductivity of self-dual heterophase systems is 
a nonuniversal function of partial parameters of phases and depends on the 
structure of their inhomogeneities. 

2. Symmetry and duality properties 

Since the inhomogeneous heterophase systems can be very different, we begin 
with a discussion of some general properties of the effective conductivity <r e 
of two-dimensional N -phase systems with the partial conductivities Cj [i = 
1,2, N) and concentrations Xi Q2 1 x, = 1). For a brevity, we will denote 
a set of all partial parameteres, for example, x«, (i = 1, ■■■,N) as {x}. It was 
already noted above that there are continuous heterophase systems and their 
discretized RRN type models. Some properties of these systems are different 
due to a discretization procedure, which is approximate and rather ambiguous. 
For this reason RRN type models are not always a good approximation for 
continuous media. 

The effective conductivity of N -phase systems must be a homogeneous (a 
degree 1) function of partial conductivities Ui and satisfies the following natural 
boundary conditions: 

1) a e of iV -phase system with some equal partial conductivities must reduce 
to er e of system with smaller number of phases and the concentrations of the 
phases with equal conductivities must add; 

2) it has not depend on partial <r, and must reduce to the effective conduc- 
tivity of (N — 1) -phase system, if the concentration X4 = ; 

3) it must reduce to partial <x; , if Xj = 1. 

As it follows from the general principles a e must satisfies also the next 
inequalities (see, for example, 

(a- 1 )- 1 <a e <{a), (5) 
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where (...) denotes the averaging with the distribution function 

N 

P(a)=J2 x i S (<r-*i)- ( 6 ) 
1 

In two-dimensional heterophase systems the effective conductivity a e satisfies 
also at arbitrary values of phase concentrations the reciprocity relation, which 
generalizes the reciprocity relation of two-phase systems [5J 1151 llb| 

a e {{a}, {x})a e ({a 2 Ja} 7 {x}) = <j 2 , (7) 

where <tq is an arbitrary constant conductivity, characterizing a reciprocity 
transformation, and a e is the effective conductivity of the reciprocal system, 
connected with the initial one. This relation is a consequence of a local duality 
between potential and divergenceless vector fields in two dimensions. The value 
co 2 /ci can be named as inversed to <j% . A structure of the reciprocal system 
depends on the type of initial system. For continuous heterophase media it 
coincides with the initial one turned on 7r/2, and for isotropic systems coincides 
with the initial one (Hj • For RRN type models the reciprocal system is the RRN 
model on the reciprocal lattice, which, in general, does not coincide with the 
initial lattice even for isotropic models. For example, in two dimensions there are 
three isotropic regular lattices: self-reciprocal square and mutually reciprocal 
triangle and honeycomb. For self-reciprocal systems the relation (7) takes the 
form 

a e ({a}, {x})a e ({a 2 /a}, {x}) = a 2 , (8) 

which puts some constraints on a possible functional form of a e . The relation 
(8) can be written also in the other form 

^({cto/ct}, {x}) = <j Q 2 /<j e ({<j}, {x}), (8a) 

which means that a e at the inversed values of the partial conductivities er^ is 
equal to the inversed a e ! For this reason (8) can be named also the inversion 
relation (IR) for the effective conductivity. The IR has a simple physical sense, 
connected with two possibilities of a writing the Ohm law in terms of conduc- 
tivities a and in terms of resistivities p = l/er : an effective conductivity of a 
self-reciprocal system as a function of partial conductivities must coincide with 
an inversed resistivity as a function of corresponding partial resistivities (i.e. 
inversed partial conductivities) pi = \jui . Thus, the DR (and IR) connects 
er e at inversed partial conductivities, but at the same partial concentrations. 
An existence of the duality relation (8) is very important, because it reflects an 
existence of some hidden symmetries, which, in principle, can help to solve the 
problem (see, for example, [T7)1. 

Another very important property of inhomogeneous heterophase systems is 
connected with their geometry and structure. For random systems with the 
symmetrical distribution function from (6) and for other heterophase systems, 
having all phases with the same properties, a e must be a symmetrical function 
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of a set of pairs of partial arguments (t7j, Xi) . For example, it fulfils for various 
symmetrical regular heterophase systems and for all RRN models, but not for 
systems with unequal fixed forms of different phases . As a symmetrical function 
er e must be invariant under all permutations 

a e ({a},{x})=a e (P ij ({a},{x})), (9) 

where Pij is the permutation of the i-th and j-th pairs of arguments (<7j, Xj) <-> 
(c7j , ) . We will call the self-reciprocal symmetrical inhomogeneous systems 
the self-dual (SD) ones. Their effective conductivity must satisfy (8) (or (8a)) 
and (9). Below we will consider mainly the SD systems. By tradition and for 
a clarity we will conserve the name "a duality relation" for (8). In the general 
form (8) it takes place for all two-dimensional continuous inhomogeneous media 
as regular as random, and for square lattice RRN model. The DR for N > i 
in general cannot be represented in a form connecting a e at adjoint concentra- 
tions, because an inversion transformation can be reduced to the permutation 
transformation only for 2-phases systems. 

An existence of the duality relation (8) contains a lot of information about 
<7 e . One can obtain from it all known general results for a e . For example, 
differentiating (8) on the parameter ao, one can show that a e satisfies the 
equation 

N 

J2^ i a e ({a},{x}) = a e ({a},{x}). (10) 
l 

It means that <r e is a homogeneous function of a degree 1 of partial conduc- 
tivities <Tj. For this reason, on the one hand, the DR for a e can be written 

as 

a e ({a},{x})a e ({l/a},{x}) = l, (86) 

from which it follows that effectively a e is a dimensionless quantity in two 
dimensions. On the other hand, it can be represented also in the form 

a e ({a},{x}) = a s f({z},{x}), (11) 

where a s is some normalizing factor, its choice is determined by a convenience of 
consideration, and z% — a/ ' a s . For our purposes it will be convenient to choose 
it in the form <r s = (sjv) 1 /^, where sjy = Y\f = i &i is the N-th elementary 
symmetric function. Then z% = ^/(sjv) 1 ^ , Ili=i z» = 1 an d a function / is a 
symmetrical function of N pairs (xi,Zi) (only N—l of them are independent). 
This choice conserves a symmetry of the dimensionless function / and works 
well for all 0$ ^ 0. Since the inversion transformation I aa acts on Zi as (see 
Appendix A6) I ao Zi = \jz, h the DR for the function / in terms of variables z 
has the following form 

f({z},{x})f({l/z},{x}) = l. (12) 

The effective conductivity of self-dual TV -phase systems, when all phases have 
the carriers with the same sign of the charge (we will confine ourselves in this 
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paper by this case), can be represented in a more symmetrical form. Introducing 
the new variables &j = ln£j, Yli 6j = 1 one obtains that a e of such systems 
is a pairwise symmetrical function of two similar sets of parameters {&} and 
{x}. 

One can also obtain from (8b) in a weak inhomogeneous case the general 
Landau-Lifshitz formula 1. . Representing the partial conductivities in this limit 
in the form <Tj = <r(l + Sai) (Yli XiSai — 0), where a — (a) is a mean value of 
cr e , and substituting them into (8b), one obtains the Landau-Lifshitz formula 
£Q, which has in a case of TV -phase systems the form 

a e ({a},{x}) = a^l-±J2 x * S(x n- ( 13 ) 

This means that the formula (13) is a consequence of general properties of self- 
dual heterophase systems. 

3. Fixed points and exact values for self-dual systems 

There is another possibility to use the DR (or IR) of inhomogeneous self- 
dual systems. Giving to uq different values, one can obtain various relations, 
connecting a e at different values of its arguments. The DR (8) can be used also 
for obtaining the exact equations for <r e ^31 • F° r this one must find such sets of 
partial parameters ({c'}, {x'}) , which are invariant under the combined action 
of some inversion transformation I aa and a set of some permutations P , i.e. 
the fixed points of all possible such combined actions. But, in contrast with a 
binary case N = 2, this cannot be done for arbitrary partial conductivities cr;, 
because, in general, a parameter <7q, for which I aa interchanges the set of all 
&i into itself for arbitrary <7j, does not exist. Such parameter can exist only 
for some special values of <7j. Then one can obtain an exact equation for er e at 
some special values of its arguments by the substitution of this combination of 
transformations in the duality or inversion relations. The exact equation will 
have a form 

a 2 e ({a'},{x'})=al (14) 

here Co is a parameter of the involved inversion transformation. The corre- 
sponding exact equation for function / has the form 

/ 2 ({z'},Lx'}) = l, (14a) 

i.e. the exact value of / at any fixed point is equal to 1. Consequently, the 
possible exact values for a e at these fixed points must satisfy the next equalities 

a e ({a'},{x'}) = a Q = (s N )^ N . (15) 

If the inversion is the interchanging transformation Jjj (A2), then 

c e = ctq = ^/o~W~ = (s N ) 1/N , i ^ j. (15a) 
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Fig.l. (a) A schematical plot of possible values of en, (i = 1,2), satisfying the relation 
(TK72 = do with some fixed value exo ; (b) a concentration phase diagram for N = 3 , 
the bisectrisses correspond to the pairwise equal phase concentrations. 

Thus, (15,15a) give a general form of the exact values of a e in these fixed points. 
The last equalities in (15) and (15a) give a strong constraint on possible fixed 
point values of C7j. 

1. 3-phase case at arbitrary concentrations. 

In order to find out possible fixed points let us consider first the case N = 3. 
Remember that in case of TV = 2 such fixed point exists for arbitrary oi (i = 
1,2), but only at the equal phase concentrations x\ = x 2 = 1/2, when the 
interchanging inversion transformation J12 can be compensated by the permu- 
tation of the pairs of arguments. This gives the known Keller - Dykhne result 
(3). The similar fixed point in case of N — 3 was found by Dykhne at equal 
phase concentrations x\ = x 2 and the additional constraint on er 3 

°3 = a = (16) 

which ensures a conservation of er 3 under the inversion. The corresponding 
exact value for a e coincides with (3). 

Now we note that the Dykhne result (16) for N — 3 admits a generalization 
on the case of equal concentrations only of two first phases x\ = x 2 = x, where 
x can be an arbitrary concentration, satisfying only the normalization condition 
2.x + £3 = 1 . In this case, due to (16), the pair of arguments (#3, 03) remains 
unchanged under the inversion transformation J12 and the permutation of two 
first pairs of arguments. Thus, the corresponding exact value for a e coincides 
again with (16). Of course, one can choose the pairs with equal concentrations by 
various ways. For example, if the phases 1 and 3 have the equal concentrations 
x and <7q = <7i<73 = of, then one obtains 

cr e = \/<y\az = gi- (17) 
Analogously, when x 2 = x 3 and ctq = 0-203 = a\, one obtains 

c e = y/wva = ci- (18) 
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Fig. 2. A schematical plot of the intersection of a plane defined by the relation Yli Oi = 
a c , where a c is some constant, with the fixed point surfaces Oi<jj = a\, i =fc j / 
k, i, j, k = 1, 2, 3, in the a -space. 

The corresponding lines of equal concentrations of different phases on the con- 
centration phase diagram are shown on Fig. lb, where they coincide with the 
bisectrisses of the triangle. Their unique intersection point corresponds to the 
case of equal concentrations of all phases. In the a -space the FP conditions 
(16)-(18) determine three surfaces. Their intersections with a plane, defined by 
the condition = a ci where a c is some constant, are shown on Fig. 2. 

Here all curves also intersect in one point, corresponding to the equal partial 
conductivities. 

It follows from above results that er e remains fixed, when Oi belong to one 
of these surfaces and the concentrations of the phases change along the cor- 
responding lines! It is worthwhile to note that now exact value of a e at all 
equal phase concentrations x\ = x<i = Xs = 1/3 is not universal, contrary to 
the N = 2 case, and depends on the constraint put on the m. The form of 
these possible exact values as a square root of product of two partial conduc- 
tivities or a value of the conductivity of the third phase is determined by the 
form of the interchanging inversion transformation and by the structure of the 
corresponding fixed points. 

2. General N -phase case. 

We show now that the DR (8) for N -phase systems admits a generalization 
of these results on the N -phase systems. 

Of course, for systems with N > 3 there are similar fixed points, when 
Xi = Xj and all other 07 (I ^ are equal between themselves. But this case 
is trivial, since it reduces to the 3-phase case. Fortunately, the new possibilities 
appear for N > 3. 

Let us consider firstly a case when N is even N — 2M. Then a new 
fixed point is possible, which corresponds to M pairs of phases with equal 
concentrations and equal products of the corresponding conductivities. For 
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example, the fixed point 

Xli-\ = X 2 i = Hi, Oq = 02,-10-2, (i = l,...,M) (19) 

is possible. But all Ui in the phase pairs must be different. In other words 
the conductivities of these pairs correspond to M different points and their 
inversed partners on the curve (see Fig. la) 

(Tq = 0102. 

The effective conductivity has in this case the next exact value 

a e = a = ^/a 2 i-i<T2i (i = l,...,M). (20) 

Note that the equal concentrations of the phase pairs y a (a = 1, ...,M) can 
be arbitrary, except the normalizing condition ^Va — 1- The exact value 

(20) does not depend on y a and ensures again the correct values for a e in the 
limits y a — > 0, 1 (a = 1,...,M). Of course, the similar fixed points exist for 
other ways of the partition of N phases into the pairs with equal concentrations. 
A number of such points is equal #2M = (2M — 1)!!. The exact values of a e 
in these points have always the same general form 

<r e = <r = y/uWj {i < j, i,j = 1, -, N), (21) 

where possible pairs (ij) correspond to the phases with equal concentrations 
in the corresponding fixed points. 

Now we consider a case when N is odd N — 2M+1 . In this case a new fixed 
point is possible if the above conditions of the even case (19) are supplemented 
by a condition on <72M+i analogous to the Dykhne condition for N — 3 case 

°2M+1 = °0 = cr 2i-lO'2i (i = l,...,M). (22) 

The corresponding effective conductivity has again the exact value (20), coin- 
ciding with the exact value of system with N = 2M . Of course, the other fixed 
points related with various ways of choice of the phase pairs with equal con- 
centrations are possible. Their number is equal to #2Af+i = (2M + 1)#2M = 
(2M + 1)!!. A general form of the exact values of o e remains the same as in 
the even case with an additional equality 

a e = <T = <Jk = ^/OiOj, i^j^k, i<j (i, j, k = 1, N). (23) 

where the pairs ij correspond to the phases with equal concentrations and k 
corresponds to the unpaired phases. 

It follows from (19)-(23) that the FP conditions determine now in the x- 
space M — 1 -dimensional hyperplanes for even N — 2M ( M -dimensional 
hyperplane for odd N = 2M + 1 ) and M -dimensional ( M + 1 -dimensional) 
hypersurfaces in the a -space. Again a e remains constant, when a, belong to 
these hypersurfaces and the concentrations X{ change along the corresponding 
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hyperplanes. One can see from (21), (23) that the exact values of a e at the 
point of all equal concentrations cc, = 1/N are also nonuniversal. 

It can be shown that only the fixed points of type (19), (22) and their vari- 
ous permutations are admissible for dual transformations of the form I aQ Y[ Pb, 
when the action of the inversion transformation is equal to (or compensated by) 
the product of permutations F = f]Pj, (see Appendix B). 

4. EM approximation for random N -phase systems 

Now, let us consider the EM approximation for the self-dual N -phase square 
lattice RRN model (we assume it also as the EM approximation for all self-dual 
N -phase systems). We will show that it has a very complicated form, and even 
to check a validity of the obtained exact results in it is not a simple task. The 
EM approximation gives the following equation for <r e El Ej 



N 

E 

i 



xt -^-^ = 0. (24) 



It can be written as an algebraic equation of the iV -th order 

N N f 

W(a e ,<Ji, ...,cr N ,xi, ...,x N ) = '^x l (a e - cr,-) JJ (a e 

i=l k=l 

r N | „ ff- 



a 1 ; + axcr^" 1 + ... + ajv-io-e + a N = 0, (25) 
dfe = s fe - 2s k = -s k + 2s k (k = 1,...,N- 1), a N = -s N . (26) 



Here 



N N i k _ 1 N N 

. — > x — , / r . ^ . — , / 



»fc 



n ^ > §k = e Xi e n ^ > ( 2? ) 



ii<i2<...<ife-i h=ii t=l h<i2< - <ik ii=ii 



where all X)'i 1 <j 2 < <i k do not contain i\ = i, and for k = N — sn, Sn — 
0. The coefficients a k are the homogeneous functions of <Ji of degree k and 
the linear functions of concentrations Xi. For this reason they are related only 
with the corresponding symmetrical functions s k and their "averages" s k or 
§k (since there are different ways to define averaged s k )■ The eq. (25) must 
have one physical solution a e > 0, which depends on all symmetrical functions 
and their " averages" . 

It is easy to check that a e , determined by the equation (25), satisfies (8a), 
since the equation (24) for inversed partial conductivities has the next property 
(for all Oi ^ ) 



N , , N 

rr ~ / rr - - rr ■ 

■r, = - > = 0. (28) 



\ ^ <7 e — <7 /(Ji __ _ V ^ gp/cTe — CT 



1 
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Moreover, since ajv = —sn one can pass directly to the similar equation for / 
by dividing all its terms on sn 



f N + a 1 f N ~ 1 + ... + a N - 1 f-l = 0, (29) 

where now the coefficients are expressed through projective variables z% 

a i = a i ({z}) = a i /(s N ) i/N , i = 1, N - 1. (30) 

It is possible to write out the physical solution of (29) in an explicit form only 
for N = 3,4, but even in 3 -phase case it has a complicated form, which doesn't 
permit to investigate it in detail or to use it for description of any experimental 
or numerical results (see Appendix C, where it is written out for a completeness). 

Here we will show that the EM approximation reproduces the exact values 
from section 3. Firstly we consider the case of equal concentrations Xi = 1/N. 
Then one can see that the coefficients and / must depend only on symmet- 
rical functions s k (k = 1, ...,N — 1) 

N — 2k 

a-k = — — s fe , / = /(si,...,Sjv-i), (31) 

Due to the transformation rules of §k (A10) the transformed / in the duality 
relation (12) has the form 

If(h, ...,sjv-i) = /(sjv-i, (32) 

Thus, one must have at the fixed point the next equalities 

s N - k = s k , k = 1, ...,N — 1. (33) 

It is easy to check that they are satisfied at the fixed points (19,22) and that 
the coefficients dk satisfy at these points the equations 

a-N-k = -Ofc- (34) 

It follows from (34) that for even N = 2M the coefficient a M = 0. Conse- 
quently, one has from (29), (34) that the polynomial in (29) can be represented 
in the factorized form 

(/-l)P JV _ 1 (/)=0. (35) 
where a polynomial Pjv-i(/) of degree N — 1 has the next form 

M i-1 M-l M-2 M-2 

ftr-i(/) = E / 2M_< E a ' + / M_1 E a ' + (E/'KE a '), N = 2M; 

»=1 i=0 i=0 1=0 1=0 

M i M-l M-l 

P J v-i(/)=E/ 2M " i E a ' + (E/')(E^), JV = 2M + 1. (36) 

i=0 (=0 (=0 (=0 
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Since all coefficients of Pjy-iCf) are positive, one can see from (36) that equa- 
tions (29) and (35) have always only one physical solution at any fixed point 
(19,22) 

f = l, a e ^{s N ) l ' N = ^] (i^j). (37) 

which coincides with the corresponding exact value from (15), (15a) or (20), (21). 

Now we pass to the case of the FP with arbitrary pairwise equal concentra- 
tions (19), (22). Firstly, we note that, due to independence of coefficients &k, 
one can write any solution / of (29) as a function of 

f = /(ai, ...,ajv-i) (38) 

Let us consider how the coefficients &k {k — 1, ...,N — 1) transform under 
inversion transformations. As it is shown in the Appendix A the transformation 
rules for at arbitrary concentrations are 

Ia a k = -a N - k , (k = 1, N- 1). (39) 

It follows from (39) that at any FP ({z'}, {x'}) one must have again 

a k = -a N „ k , (k — 1, ...,N — 1), (40) 

what reduces the polynomial in equation (29) to the factorizcd form 

(f-l)P N - 1 (f,{x'})=0, (41) 

where now the polynomial Pn-i(J, {%'}) depends on concentrations {x'}. Due 
to this factorization, eqs.(25) have at the FP (19), (22) the required solutions of 
the form (21), (23), which are the reduction of the corresponding exact values 
(15) at these points. 

Thus, we see that in the EM approximation er e of self-dual iV -phase sys- 
tem is a solution of the A-th order equation, which is an algebraic, homo- 
geneous, a degree 1, function of the set of 2N — 1 independent arguments 
{<r} = (ui, on) and {x} = (x\, xn), ( x i — 1 ) symmetric under the 
permutation of pairs of arguments (ai,Xi) and simultaneously satisfying the 
inversion relation (8a) and the exact values (23). The solving of the equation 
(25) in general N -phase case is an interesting mathematical task due to its 
relation with the symmetrical and inversion groups. Unfortunately, it is not so 
useful for the understanding of physics of the problem, since already in 3-phase 
case the EM approximation gives for a e a complicated algebraic expression, 
which is cumbersome and practically untractable (see Appendix C). This can 
be explained by the fact that this solution depends on all symmetrical functions 
Si, (i = 1, 2, N). From the other side, it was shown in the binary case N = 2, 
that there are other simple solutions of the DR, which, though differ from the 
EM approximation expression, but also satisfy all necessary requierements and 
can describe a e of some self-dual continuous inhomogeneous systems for not too 
large inhomogeneities with the same accuracy as the EM approximation jll| . 
In the next sections we will show that analogous solutions and systems exist for 
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arbitrary N and that they have the simple forms, convenient for approximate 
description of the effective conductivity in a wide range of parameters. Their 
simplicity can be explained by the next circumstence: they are related only with 
one symmetrical function and its averaged versions. 

5. Solutions of inversion relation 

1. General solution. 

In this section we consider the possible functional forms compatible with 
the duality (8) and/or inversion (8a) relations. Due to its generality, the DR 
cannot define the effective conductivity unambiguously and completely. The 
diversity of the possible functional forms of a e can reflect a variety of all possible 
inhomogeneity structures of random inhomogeneous systems. For this reason 
it is very desirable to find all possible functional forms compatible with the 
duality relation or to propose a general method of their construction. Since 
(8) and (8a) take a place for all self-reciprocal systems, the functional forms 
compatible with the DR will be applicable for all self-reciprocal systems. The 
different self-reciprocal systems are distinguished by corresponding additional 
functional properties. The self-dual case corresponds to the additional pairwise 
permutation symmetry of <r e . The duality and/or inversion relations (8), (8a) 
must be supplemented by the boundary conditions 1) - 3) enumerated above in 
the Section 2, which turn out rather restrictive. 

For further notation simplification it will be convenient to consider the DR 
in the form (8b). Then, taking into account the properties of the inversion 
transformation (A3), one can write out formally a general solution of (8b) in 
the following form 

a e ({a}, {x}) = g{{a], {x})/g({l/a}, {x}), (42) 

where g is an arbitrary function. But, one needs to remember that a e must be 
a homogeneous function of a degree 1 and obey additional requirements. The 
simplest way to satisfy the homogeneity condition is to require g to be also a 
homogeneous function. Then, if g is homogeneous of a degree k, one can show 
that cr e is 

<7e(M, W) = (ff(M. {x})/g({l/v}, {x})) 1,2k . (43) 

In order to obtain some a e for self-dual systems a function g must be a pairwise 
symmetrical function and satisfies the boundary conditions l)-3). 

It is useful to consider a general solution in the projective variables {z}. 
Then for a function / from (11) one obtains from (43) the general solution 

/({*}, {x}) = (g({z}, {x})/g({l/z}, {x})) 1/2k , (44) 

where a function g is a pairwise symmetrical and is homogeneous function of 
a degree k of {z}. 

2. Some basic solutions. 
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Now we consider some simple solutions of the inversion relation for self-dual 
systems. 

1) As a simplest admissible pairwise symmetrical function one can take the 
mean conductivity g({a},{x}) = (cr) = si = J2 i=1 Xi(Ji, which is a pairwise 
symmetrical homogeneous function of a degree 1. It satisfies also the boundary 
conditions. The corresponding function of the inversed partial conductivities 
g({l/a}, {x}) = (cr -1 ) = J2i x i/ a i — sjv-i/sjv- Then, substituting these func- 
tions in (46) one obtains 



^(M.M) -VM^ 1 }- 1 = (s 1 s N /s N . 1 ) 1/2 = (s N )^ N (Ii/Ijv-i) 



1/2 



_ [ Li=i^ilifc=i_gj \ a ,x 

\ Ei=l X i II k=l °k ) 

It is clear from the form of the solution (47) that being a geometric average of 
(cr) and (cr -1 ) -1 it satisfies the inequalities (5). It is easy to check that (47) 
satisfies the exact values at the FP. 

2) Analogously, one can take as g the averaged k -th moments of a conduc- 
tivity (cr fe ) (fc = 2,3,-)> 

TV 

gk({a},{x}) = (a k )=J2^l (46) 
l 

which are the pairwise symmetrical homogeneous functions of a degree k, sat- 
isfying boundary conditions. They can be considered as the "averaged" version 
of the newtonian symmetrical functions Nk = of. The corresponding 
functions of the inversed partial conductivities are 

N N N N 

3fc ({i/cr}, { X }) = £ Xi /o* =5>n' n ^ = (47) 
i ill 

where = Ef^Il'f CT f- Substituting (46) and (47) into (43) one obtains 

the solutions 

<#>({*},{*}) = (^^/^li) 172 " - (^) VAr • (48) 

These simple algebraic solutions can be considered as the basic solutions, cor- 
responding to the k -th moments of a conductivity. They satisfy all necessary 
requirements. It is worth to note here that the higher moments of conductivity 
fluctuations are already accounted, due to the DR, even in the simplest solution 
(45). 

In order to see the importance of the boundary condition requirements one 
can consider as the functions gk the "averaged" functions Sk or §k (k = 
2, 3, N — 1), appearing in the EM approximation in the section 4 and defined 
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by (27). They are homogeneous, pairwise symmetrical and, in some sense (see 
(A7)), inversed to each other. Substituting these functions into (43) one obtains 
the following solutions 

4 k) (Wh W) = (s k s N /§N-k) 1/2k = 

(sn) 1/N (s k /sN-k) , (k = 2,...,N-l) (49) 

or 

4 fe) (M,{*}) = {~s k s N /s N - k ) 1/2k = 

/~ _ \ I /2k 

{s N ) 1/N [h/sN-k] , (k = l,...,N-2). (50) 

Unfortunately, these solutions do not satisfy the reduction and boundary con- 
ditions from the Section 2, because, for example, for a% = U2 

(jV) (N-l) (N-l) ~(N) (N-l) (N-l) 

S k ^ S k +°lS k _i , S k -> S k +(T 1 S k _ 1 , 

where the upperscript (N) denotes the corresponding N -phase system. It 
follows from these formulas that under the reduction o\ = &i the functions 
s k and § k (except si and §n-i ) relate with the corresponding lower 
functions. For this reason the solutions (51) and (52) are not physical ones. 
Since the "pure" symmetrical functions reduce by the same law, the coefficients 
functions a k from the section 4 also reduce for ax = &2 by the same law, and 
because the solutions connected with them are also nonphysical ones. 

3) It is not possible to define in this way "an average" for N -th symmet- 
rical function. Fortunately, there is a simple pairwise symmetrical exponential 
function, generalizing the N -th symmetrical function, the monomial fL =1 cr^ ■ 
It can be represented as an exponential of the simplest pairwise symmetrical 
function 

N 

g = exp^ biXi, h = lntTi. 

i=l 

It is easy to see that it is invariant under the joint inversion of the partial 
conductivities and of the function itself 

g- l ({l/a},{x}) = g({a},{x}). 

For this reason its substitution into (43) reproduces itself 

N 

^ iv) (w,w)=n< i - ( 5i ) 

The solutions (45) and (51) generalize two solutions found earlier for 2-phase 
systems Since they satisfy the inversion relation, they automatically have 

all necessary properties determined by this relation. For example, it is easy to 
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check that (51) reproduces at the corresponding fixed points the exact values 
from (20), (21), (23) due to the normalization condition for concentrations 



°* = IT^ = °* = (sn) 1/N , (N = 2M), (52) 



i=l 



M 

<J, = o 1 - 2 ^" y > l[ a 2 v ' = a = (sn) 1/N , (N = 2M + 1). (53) 

i=l 



It is interesting that in this case a e at the point of all equal concentrations 
Xi = l/N always coincides with s]^ N , (i.e. for arbitrary o"j). The solution 
(51) satisfies the inequalities (5) due to the Holder inequalities 

N N N N 



n °i* ^ n °i xi ^ J2 ^ 



N 

E 



x l = 1, Oi> 0. (54) 



l 



The physical models with effective conductivities equal to the simplest solu- 
tions (45) and (51) will be represented in the two next sections, where the other 
physical properties of these solutions will be also considered. The basic solutions 
(48) can be used also as constructive blocks for more sophisticated solutions. 
For example, let fi (i — 1, 2, ...n) are the solutions, satisfying all requirements. 

Then / = (J|™ /i) 1 ^™ is also the solution, satisfying all requirements. The phys- 
ical models related with the higher moments solutions (48) remain unknown at 
this time. One can also use as constructive blocks the solutions (49), (50), 
but this is a complicate task to arrange them in a combination, satisfying all 
necessary requirements. The EM approximation case confirms this. 

6. FMSA approximation and functional equation 

In this section we show that the solution (51) can describe a e of continuous 
inhomogeneous self-dual systems with compact inclusions of different phases. A 
schematical picture of such systems is shown on Fig. 3a. In order to do this we 
will generalize an approximate functional equation for er e of two-phase systems 
from on A -phase systems and solve it. 

To obtain this functional equation we can use a method called the FMSA 
approximation and proposed in |12| . One can see that all arguments represented 
there in favor of the FMSA approximation hold also for random inhomogeneous 
heterophase systems with any finite number of phases N > 3 . 

To do this let us consider the auxiliary square lattice with the squares of a size 
I ^ Imax j where l max is a higher boundary value of l m (x) , the maximal size of 
inhomogeneities of the system (which can depend on the phase concentrations) , 
i.e. l m (x) always remains smaller than l max . It is assumed that the initial 



16 




1 


M 


1 


1 


2 


2 


1 


1 


2 


2 


2 


l 


2 


2 


2 


2 


2 


1 


1 


2 


1 


2 


1 


2 


1 


2 


1 


1 


2 


1 


2 


1 


2 


2 


1 


2 


2 


1 


2 


2 


1 


2 


1 


1 


1 


2 


1 


1 


1 


1 



(b) 

Fig.3. a) A schematic picture of the elementary square with a size I ;§> l m of one phase 
with compact inclusions of the other phases, b) the auxiliary lattice of the model, the 
numbers 1,2 denotes squares with the corresponding sets of concentrations {a;}i,2 • 



inhomogeneous system has such l max , at least, it exists in some range of the 
partial concentrations. We require that the squares of the auxiliary lattice have 
the effective conductivities corresponding to two different sets of the concentra- 
tions {x}i and {x} 2 with equal probabilities p = 1/2 (see Fig. 3b). Then, 
implementing to this lattice the exact Keller-Dykhne formula (3) and supposing 
for a self-consistency that a e of this lattice must be equal to the same effective 
conductivity at the averaged concentrations {x} — {{x}i + {x} 2 )/2, one can 
obtain the approximate functional equation for the effective conductivity of this 
system as a function of all partial concentrations {x} 

a e ({a}, {x}) = V ^e({o}, {x}i)<T e ({a}, {x} 2 ), 

{x} = ({x} 1 + {x} 2 )/2, (a^O). (55) 

The eq.(55), in contrast with the DR, connects a e at different concentrations, 
but at the same conductivities. Strictly speaking, the eq. (55) takes place only 
for the ranges of Xi, where l max exists. These ranges can be enlarged by 
taking the larger I. The boundary conditions for the equation (55) follow from 
the conditions 1) - 3) of the Section 2 

a e ({cr},Xi = 1) = <7j, a e ({a},Xi = 1) = a t , (55') 

Let us suppose for a moment that the equation (55) is applicable for all phase 
concentrations. Then a solution of (55) has an exponential form, generalizing 
the solution of two-phase systems from El 

N N 

a e ({a},{x}) = l[a i x % = 1. (56) 

i=l i=l 

It coincides exactly with the solution (45) and has the simplest form compatible 
with both symmetry groups (symmetrical and inversion ones) of the problem. 
We see that constrainting himself by the RRN type models or by the EM approx- 
imation and the corresponding algebraic equation one cannot obtain the more 
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simpler exponential solution. At equal concentrations of all phases x% = 1 /N 
one obtains from (56) 



N 

a e ({a}, {x}) = J] ^ 1/N = (sn) 1/N - (57) 

i=l 

This formula generalizes the exact Keller ~ Dykhne formula (3) for systems with 
N = 2 on the case of random inhomogeneous systems with arbitrary N phases 
and the exact values (23) for N -phase systems at the FP on all Cj ^ 0. 

But, as it was mentioned in the section 3, the formulas for a e at Xi — 1/N 
are not universal. For example, it differs from the corresponding formula in 
the EM approximation (see Appendix C). Below we will obtain another formula 
for o~ e of the random N -phase model with other structure of inhomogeneities, 
which differs at Xi — 1/N from both these formulas. 

It is easy to check that in the weak inhomogeneous limit, when all partial 
conductivities ai can be represented in the form 

(Xi = cr (l + 8<Ti), 

where Co is a constant homogeneous conductivity, the formula (56) gives up to 
the second order on 5a,; 



°~e = ^0 JJ(1 + 8<7i) X ' W (T 



N 1 



which reproduces the general Landau-Lifshitz formula (13) for ao = (c), when 
^XiSdi = 0. The formula (56) means that lner is a self-averaging quantity in 
this approximation, because it can be represented in the form 

N 

a e ({x}, {<t}) = exp(y~^ Xj In aj) = exp(lncr) = exp(b). (58) 

i=l 

This property noted firstly by Dykhne for the two-phase systems at equal con- 
centrations x = 1/2 JTTJ was established later in the weak localization theory 
of conductivity of two-dimensional disordered systems ^2] . Since the solution 
(56) reproduces the exact values at all FP, it can be used in a more wide range 
of concentrations than it was assumed under the deducing of the eq.(55). We 
expect that it works well for Cj > at all concentrations for all self-dual in- 
homogeneous systems, to which the FMSA approximation is applicable. For 
example, it can be applyed to the inhomogeneous disordered systems where the 
regions with inhomogeneously distributed impurities can be described effectively 
as the regions of different phases. One needs to note that the formula (53) for 
arbitrary partial concentrations have been proposed earlier as a result of the 
extrapolation of the analogous formula for o~ e of some special models of ran- 
dom isotropic polycrystal systems at discrete set of concentrations on arbitrary 
concentrations [T5). 
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Fig. 4. (a) An elementary square of the model with a vertical orientation of the layers 
of different phases; (b) a lattice of the model, the small lines on the squares denote 
their orientations. 



7. Random parquet model of heterophase systems 

Now we will construct a hierarchical model of two-dimensional isotropic ran- 
domly inhomogeneous heterophase system, using the composite method intro- 
duced in |21 112j , and show that its effective conductivity a e obtained in the 
mean field like approximation for arbitrary values of the phase concentrations 
Xi coincides with the solution (45) of the DR. 

Let us consider the following two-dimensional model. There is a simple 
square lattice with the squares consisting of a random layered(or striped) mix- 
ture of N conducting phases with constant conductivities <7j,i = 1, ...,7V and 
the corresponding concentrations Xj. A schematic picture of such square is given 
in Fig. 4. 

The layered structure of the squares means that the squares have some pre- 
ferred direction, for example, along the layers. Let us suppose that the directions 
of different squares are randomly oriented (parallely or perpendicularly) rela- 
tively to the external electric field, which is directed along x axis. In order 
for system to be isotropic the probabilities of the parallel and perpendicular 
orientations of squares must be equal or (what is the same) the concentrations 
of the squares with different orientations must be equal pu — p± = 1/2. 

Such lattice can model a random system consisting from mixed phase regions, 
which can be roughly represented on the small macroscopic scales as randomly 
distributed plots with the effective "parallel" and "serial" connections of the 
layered or striped N -phase mixture (Fig. 4). The lines on the squares denote 
their orientations. This structure can appear, for example, on the intermediate 
scales when a random medium is formed as a result of the stirring of the N - 
phase mixture or as a result of spontaneous generation of randomly distributed 
regions of N striped phases with different orientations of stripes near the phase 
transition point. For external electric field directed along one of the preferred 
orientation the corresponding averaged parallel and perpendicular conductivities 



19 



of squares m | (a;) and a± (x) are defined by the following formulas 

N 

i 

a x ({a},{x}) = ^Xi/a^j = (a" 1 )- 1 . (60) 

Thus we have obtained the hierarchical representation of random medium (in 
this case a two-level one). On the first level it consists from some regions (the 
squares with two possible stripe orientations) of the random mixture of the 
N layered conducting phases with different conductivities 0$ and arbitrary 
concentration. On the second level this medium is represented as a random 
parquet constructed from two such squares with different conductivities m | and 
a± , depending nontrivially on concentrations of the initial conducting phases, 
and randomly distributed with the same probabilities p, = 1/2 (Fig. 4). This 
representation allows us to divide the averaging process into two steps. Firstly, 
we average over each square and obtain for them crn or <jj_. Then, on the 
second step, the averaging is done over the lattice of squares and the exact 
Keller - Dykhne formula (3) is used. As a result one obtains for the effective 
conductivity of the N-phase random parquet model the following formula, which 
is applicable for arbitrary concentration 

a e ({a},{x}) = V^I=v / W(0- 1 = l tt= N l(Tk ^ lXiai ■ (6 1 ) 

It coincides with the solution (45) from the section 5 and generalizes the cor- 
responding solution for two-phase random parquet model from [OJ I12| . It is 
worth to note here that such formula for er e has been obtained earlier for inho- 
mogeneous media with a smooth log -symmetrical distribution function, when 
it is the exact one [TT]. Since such inhomogencous media correspond to the 
limit N — > oo and a smooth log -symmetrical distribution function of the ran- 
dom parquet model, the solution (61) becomes exact in this limit. Analogous 
expression for a e has been also proposed for one model of random isotropic 
polycrystal systems with a structure of inhomogeneities similar to those of the 
random parquet model |19| . 

One can also check that (61) reproduces the FP values (15), (21) and (23). 
In a case of the equal phase concentrations Xi = 1/N one obtains now 

a e ({a}, {1/JV}) = y/siSN/sN-L (62) 

Thus we see that the expressions (45), (51) and the EM approximation (Ap- 
pendix C) for the effective conductivity have different functional forms. This 
fact means that the effective conductivity of self-dual random inhomogeneous 
systems is a nonuniversal function. Moreover, it follows from these expressions 
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that the values of a e at equal phase concentrations Xi = 1/N and arbitrary 
(Tj are also nonuniversal (unlike the case N = 2). This is a consequence of 
the nonuniversality of the effective conductivity and a fact that this point at 
arbitrary partial conductivities is not the fixed point. A nonuniversality of a e 
at Xi — 1/3 for some 3- phase self-dual regular inhomogeneous systems was also 
noted earlier on a basis of numerical calculations in the paper |16| . 

One can easily find from the explicit formulas (45) and (51) the correspond- 
ing low concentration asymptotics of <r e and see that their coefficients diverge 
at ex, — > 0. This shows that the limits er,; — > can have some singularities as 
in the two-phase case [91 IT2]. 

9. Conclusions 

Thus, using the duality relation and some additional assumptions, we found 
out a general form of the effective conductivity a e of two-dimensional self-dual 
heterophase systems and written out some explicit approximate expressions for 
a e at arbitrary values of partial concentrations and conductivities. . Though 
these formulas have the different functional forms, they 

(1) satisfy all necessary requirements including symmetries and inequalities, 

(2) reproduce the general Landau - Lifshitz formula for a e in the weakly 
inhomogeneous case, 

(3) reproduce all exact values for <j e at the fixed points. 

One of them is also a solution of the approximate functional equation, and 
can describe tr e of various self-dual (or close to them) heterophase systems 
with a compact inclusions of phases or with a finite maximal scale of inho- 
mogeneities l max . We have constructed also a random parquet model of the 
random inhomogeneous medium of the layred type and have shown that its ef- 
fective conductivity in some mean field type approximation coincides with the 
solution connected with a mean conductivity. All these results show also that, 
in general, tr e of the inhomogeneous heterophase self-dual systems may be a 
nonuniversal function and can depend on some details of the structure of the 
randomly inhomogeneous regions. The obtained explicit formulas for a e can 
be used for an approximate description of the effective conductivity of some real 
random and regular inhomogeneous systems. 

Appendix A. Inversion, symmetrical functions and fixed points 

In this appendix we consider the inversion group and its action on various 
symmetrical functions. The inversion group with a parameter <tq, acting on 
one variable function f(a) , contains two elements 

I„ f(a) = f(a 2 /a), 1^=1. (Al) 

Since the effective conductivity of the N -phase system depends on N partial 
conductivities, one must consider the action of the direct products of the partial 
inversion transformations I^o ■ a% —> &o/ci, (i = 1,2..., N). Each partial inver- 
sion transformation acts only on the partial conductivity of the corresponding 
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phase and is characterized by an arbitrary parameter oq. Formally, the inver- 
sion transformation is well defined only for all <7j ^ 0. One can extrapolate 
it on the case, when some <7j = 0. In this case it transforms a system with 
<7j = into a system with a% = oo. Note also that the inversion transformation 
with parameter <r has only one fixed point <7j = <t - The i-th and j-th partial 
inversion transformations iff (iff) , having (Tq = a%aj, are very important, 
because they change oi into aj and vice versa 

J ff a i = P ij = ViOj/oi = <>j, iffaj = aiUj/uj = Gi. (A2) 

The physical sence has only the product of all partial inversion transformations 
with the same parameter a 

N 

I„ =l[l%, Il =l- (A3) 

i=l 

Since in our problem the inversion transformation acts on a space of symmet- 
rical homogeneous functions it is useful to consider its action on the basic el- 
ementary (and newtonian) symmetrical functions Sk({cr}) (k = 1, N, {a} = 

((7i, <7 N )) 

N N N i fc 

i=l i<j ii<i2<---<ik il 

h Sk = <rl k s N ^ k /s N . (AA) 

It follows from (A4) that only s n transforms into itself. Under consideration of 
the inversion transformation action on homogeneous functions it is convenient 
to pass (for cr, ^ ) to the projective variables Zi and introduce projective 
basic symmetrical functions §k 

Sk = s k ({z}) = s k /(s N ) k / N , k = l,..,N-l. (A5) 

Then we have from (A4) and (A5) 

I ao Zi = l/zu I aa s k = SAr_fe, (A6) 

i.e. the inversion transformation I ao for any cto inverses Zi and interchanges 
adjoint s^! This property will strongly simplify a search of fixed points needed 
for obtaining the exact equations for a e . Now we consider how the coefficients 
a/c (k = 1, ...,N — 1) from the EMA equation (25) transform under inversion 
transformations. Since a k = s k — 2s k = —s k + 2s kl one needs to find the 
transformation rules for s k and/or Sk- Using the definitions of s k and s k 
from (29) and one can show that 

I<r Sk = s N - k /s N , I ao s k = s N - k /s N , (A7) 
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In the projective variables Zi these transformations simplify 

I<j a Sk — SN-k, Icr Sk = §N-k, (A8) 

These relations generalize the relations for pure symmetrical functions (A6) on 
a case of arbitrary Xi. It follows from (A7) and (A8) that the transformation 
rules for a k at arbitrary concentrations are 

I<j a a k = -a N -k, (k = 1, N - 1). (A9) 
The fixed points {a 1 } of the inversion transformation I ao are defined as 

I*oW'} = W} (AW) 

Thus, at any FP one must have 

a k ({a'},{x'}) = -a N _ k ({a'},{x'}), (k = 1, JV - 1). (All) 



Appendix B. Evidence of absence of other fixed points 

One can show that only the fixed points of type (19,22) and their various 
permutations are admissible for all dual transformations of the form T(I, P) = 
[\I a Y\Pb, where, for clarity, both products are assumed to be finite, and the 
action of the product of inversions I = J| I a must be equal to (or compensated 
by) the product of permutations P = FJ P5 (we will not consider here possible, 
more sophisticated, forms of T). To do this it is useful to pass from partial 
conductivities <7j to their logarithms b, = lnc^, (—00 < bi < 00). Then the 
partial inversion transformation l£j acts on bi as a reflection bi — > B — bi, B = 
2 In (To- The inversion I aa acts on the set of variables {6} in the following way 

bi^B-bt, (i = l,...,JV) (BI) 

The action of the product of L different inversions I = nf=i gives anal- 
ogous result 

h^B L + (-l) L b t , (i = l,...,N) (B2) 

with Bl = J2s=i(~^) s B s , B s = 21ncTQ S ' ) . Since all parameters of inversion 
transformations collect into one constant Bl, it will be enough to consider 
only L = 1 or L = 2. In order to find any fixed point the result () must be 
equivalent to some permutation P of the set {b} 

{b} P = (bi p ,...,b Np ). 

Let us consider firstly a case of even L = 2. Then one obtains from (25) 

b lp =B 2 + b u (i = l,...,N) (B3) 
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The set of equations (B3) is very restrictive. If any bi remains unchanged, than 
B 2 = and all other bi also remain unchanged. It means that in this case 
1=1. If () gives a permutation of any pair bi <-> bk, (i ^ k), than it follows 
from () that again B 2 = 0. As is known any permutation P decomposes into 
product of independent chains of simplest permutations. Then, considering all 
possible chains similar to the abovcmentioncd simplest case of pair permutation, 
one can see that equations (26) take place only when B 2 = 0. In case of odd 
L = 1 equations (25) take the form 

h P =B 1 -b l , {i = l,...,N) (54) 

If (27) gives a permutation of some pair bi <-> bk, (i ^ k), then one obtains for 
bi, bk only one equality (since in this case both equations coincide) 

b l + b k =B 1 . (55) 

In case of any larger chains of permutation one obtains additional equalities, 
which have only trivial solution. For example, in case of 3-chain ( 1 — > 2 — > 3 — > 
1 ) they are 

&i + b 2 = b 2 + 63 =h+ 63 = B 1 . (56) 

It follows from (29) that this is possible only if bi = b 2 = 63 = b, B\ = 26, 
but this corresponds to the one phase case. Consequently, the equations (27) 
can take place if and only if the permutation P decomposes into product of 
pair permutations (for N = 2M ) or into product of pair permutations and one 
identical transformation (for N = 2M + 1 ) . Then one obtains in case of pair 
permutations of type &21-1 b 2i the next equalities 

b 2 i-i+b 2i = B 1 (i = l,...,M), (JV = 2M), 

b 2l -i+b 2l = 2b 2M+1 = B 1 (* = 1,...,M), (N = 2M + 1). (B7) 

All other solutions can be obtain from (30) by various permutations. These 
solutions give exactly the fixed points (19,22) after returning to the variables 
{a}. 

Appendix C. EM approximation for the 3 -phase RRN model 

For 3-phase case the solution of the EMA equation can be presented in an 
explicit form, but we will see that it is very complicated and inconvenient for 
a description of experimental data and theoretical analysis. According to the 
general results of the section 4 the equation for effective conductivity of the 
random 3 -phase model is 

a\ + dial + a 2 a e + a 3 = 0, 

ai = si - 2si = —81 + 2si, a 2 = s 2 - 2s 2 = -s 2 + 2s 2 , c = —S3. (CI) 

Here s\ = a = Y^i x i a i IS a m e an conductivity, s\ = x\(o 2 + 0-3) + x 2 {o\ + 
0-3) + x 3 (ai + <j 2 ), s 2 = xi(ai 2 + (7i3 ) + x 2 (a 12 + a 2a ) + x 3 (a 13 + a 23 ), s 2 = 
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a;i<723 + X2<J\3 + X3P12, (pik = crjCTfc). The reduced form of the EMA equation 
is 

y S +py + q = 0, (C2) 

where y = a e + ai/3, p = — af/3 + a 2 , </ = 2(ai/3) 3 — aia 2 /3 + a 3 . Three 
solutions of the equation (C2) are 

yi=A + B, 2/2,3 = -(A + B)/2 ± i(A - B)/2>/3, 

A=(-| + VQ) 1/3 , B^(-f-^) 1/3 , ^(f) 3 + (f) 2 . (C3) 

In order to find the physical solution let us consider firstly the homogeneous limit 
£71=0-2 = <73 = cr. Then, using the normalization condition for concentrations 
Si=i x i = 1 an d the degenerated values of Si ( s\ = 3cr, s 2 = 3cr 2 , s 3 = ct 3 ) 
one can see that 

a 1= a 7 a 2 = -a 2 , p=-^a 2 , q = -2(^fa 3 Q = 0, A = B = ^a, (C4) 
and consequently 

4 2 

2/1 = gO", J/2,3 = -gCT, CT e i = (7, (7 e2i3 = -(7. (C5) 

Thus, the solution j/i corresponds to the physical one. The expression (C3) 
has a complicated form in general case. Let us consider the physical solution at 
some special values of parameters, when it takes more simple, tractable form. 
1) Equal phase concentrations, Xi = x = 1/3. 

In this case all coefficients can be expressed only through Sj. Then one 
obtains for y\ 

a = Si/3, b = -S2/3, p = {s 2 + 9s 2 )/27, q = (2s 3 - 27sis 2 - 27 2 s 3 )/(27) 2 , 

Q = (729s 2 . - 4s 3 s? - 54s 3 s 2Sl - s 2 s^/3 - 4s 2 ! )/427 2 . (C7) 

Even in this case it is not so simple to check that a e satisfies the inversion 
relation. The equality of the numerical coefficients before a\ and a 2 is a 
consequence of the duality relation. The simple result for a e , coinciding with 
the exact values (16-18), can be obtained only in the corresponding FP a%aj — 
\tijk\c 2 , Xi — Xj, (i, j, k = 1, 2, 3), where djk is the fully antisymmetric tensor 
and there is no summation on (i, j, k). 2) To see this let us consider the FP with 
arbitrary equal concentrations of two first phases: <7i<r 2 =03, x\ = x 2 = x. In 
this case 

Sl = (71 + <J 2 + s/<fl<T^, S 2 = y/<Tl<J 2 Sl, s 3 = (cict 2 ) 3/2 , 

s 1 = a = x(ai + c7 2 ) + (1 - 2x)^(7ic 2 , s 2 = sfcno^o, 
a\ = S\ — 2a a 2 — y/<Ji<j 2 (—Si + 2a) = —a ly /a 1 a 2 . 
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Thus one see that the coefficients of equation (CI) satisfy the fixed point rela- 
tions (42) for arbitrary x. It means that eq.(Cl) has the solution / = 1, a e = 
{sz} 1 ^ = = 03. The similar behaviour takes place for other fixed points. 
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